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Abstract;  An  X-ray  diffraction  investigation  was  carried  out  in  an  attempt  tc 
determine  in  detail  the  atomic  positions  and  temperature  motions  in 
the  crystal  structure  of  the  tetragonal  ferroelectric  form  of  barium 
titanate.  The  basic  experimental  data  used  in  the  analysis  con- 
sisted of  the  integrated  intensities  of  99  Bragg  reflections  of  the  type 
(hkO).  The  intensities  were  carefully  measured  by  means  of  a scan- 
ning technique  using  a Geiger  counter  assembly  especially  designed 
for  this  study,  a small  (0.1  mm)  untwinned  single  crystal,  and  MoKa 
radiation.  The  intensities  were  corrected  for  absorption  and  second- 
ary extinction  effects  determined  by  special  techniques  worked  out  for 
this  particular  crystal. 

The  search  for  a structure  which  will  adequately  account  for 
these  data  was  carried  out  by  the  methods  of  trial  and  error  and 
least-squares  analysis.  A model  was  finally  defined  in  terms  of  3 
structure  and  9 anisotropic  temperature  parameters.  It  was  found 
that  in  the  general  intensity  function  for  all  reflections  there  is  an 
interaction  between  the  structure  parameter  of  an  atom  (displacement 
along  the  tetragonal  axis)  and  the  temperature  motion  of  the  atom  in 
this  direction.  Therefore,  in  order  to  determine  the  value  of  one  from 
the  observed  intensity  data,  it  is  necessary  to  assume  a fixed  value  for 
the  other.  Consequently,  within  a rather  wide  range,  it  is  not  possible 
to  determine  uniquely  the  structure  of  tetragonal  barium  titanate. 

Certain  temperature  parameters  do  not  interact  with  other  para- 
meters, and  these  are  uniquely  determined.  From  these,  reasonable 

*A  considerable  portion  of  the  work  described  here  has  been  reported 
in  Technical  Reports  Nos.  39  and  54,  Philips  Laboratories,  Inc. 

**Present  address;  U.  S.  Geological  Survey,  Washington,  D.'C. 
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values  have  been  assigned  to  the  other  temperature  motions,  and  the 
structure  consistent  with  these  assumptions  determined.  The  reli- 
ability factor  R for  this  structure  is  0.0367.  These  results  are 
compared  with  those  published  by  other  workers. 

1.  Introduction 

Barium  titanate  (BaTiO^)  is  one  of  a large  family  of  compounds  having 
the  perovskite  (CaTiOj)  structure,  or  a close  approximation  to  it.  The  ideal 
cubic  structure,  realized  by  BaTiO^  above  120°C,  is  very  simple  and  is  il- 
lustrated in  Fig.  lb.  Most  compounds  of  the  perovskite  type  actually  have 
Slight  distortions  of  this  ideal, highly  symmetrical  structure,  and  barium  titanate, 
in  particular,  becomes  polar  below  120°C.  At  this  point  BaTiOj  shows  most  un- 
usual dielectric  properties  (dielectric  constants  as  high  as  20,000  have  been 
measured),  and  becomes  piezo-  and  ferroelectric.  Because  of  its  simple 
structure,  and  convenient  chemical  and  unusual  physical  properties,  barium 
titanate  has  been  subjected  in  the  last  eight  years  to  intensive  theoretical  and 
experimental  investigation  in  many  laboratories  throughout  the  world  as  an  out- 
standing problem  in  solid  state  physics 

It  has  been  commonly  stated  that  the  polar  properties  arise  from  a 
slight  displacement  of  the  titanium  atom  from  its  symmetrical  position  in 
the  center  of  any  oxygen  octahedron  toward  one  of  the  oxygen  atoms  (in  the 
direction  of  a cube  axis).  Actually,  all  atoms  may  be  displaced,  and  a theo- 
retical understanding  of  the  properties  of  BaTiOj  must  depend  on  an  exact  know- 
ledge of  these  displacements.  At  first  glance  it  appears  that  the  methods  of 
X-ray  diffraction  should  readily  yield  this  information,  but  many  unexpected  dif- 
ficulties and  special  features  of  the  diffraction  problem  have  delayed  the  appear - 

1)  A.  R.  von  Hippel,  "Ferroelectricity,  Domain  Structure  and  Phase  Transitions 
in  Barium  Titanate,  " Technical  Report  No.  23,  Laboratory  for  Insulation 
Research,  March,  1950. 


-4- 


ance  of,  and  qualified,  the  desired  results. 


2.  Phases  and  Structure  of  Barium  Titanate 
Two  polymorphs  of  barium  titanate  are  well  known,  one  hexagonal  and 
one  cubic  or  pseudocubic.  According  to  Burbank  and  Evans^  these  struc- 
tures are  both  based  on  close -packed  barium  and  oxygen  ions,  which  are 
practically  the  same  size  but  having  different  stacking  sequences  (Figs,  la  and 
lb).  The  titanium  ions  are  found  in  octahedral  oxygen  interstices  between 
close-packed  layers.  In  the  cubic  form,  the  close-packing  is  cubic,  and  the 
TiO^  octahedra  are  linked  by  corners  in  straight  rows  throughout  the  structure, 
unlike  the  hexagonal  form  where  the  titanium  distribution  is  more  complex. 

The  cubic  form  therefore  presents  the  possibility  of  a uniform  interaction  be- 
tween Ti-0  dipoles  through  the  structure,  producing  various  polar  deformations 
of  the  cubic  structure  at  various  temperatures. 

In  the  perovskite  type  of  barium  ti+anate,  four  phases  are  known,  well 

3 4)  5) 

established  by  optical  and  X-ray  measurements,  ’ ' and  are  defined  ' as  fol- 
lows: 

I.  Cubic}  space  group  = Pm3mj  aQ  = 3.9860.& 

Cell  contents:  BaTiO^ 

Ba  in  (a),  0,  0,  0 

Ti  in  (b),  l/2,  l/2,  l/2. 

3 O in  (c),  1/2,  1/2,  0}  l/2,  0,  l/2{  0,  l/2,  l/2 
No  parameters;  stable  above  120°C. 

2)  R.  D.  Burbank  and  H.  T . Evans,  Jr.,  Acta  Cryst.  1,  330(1948). 

3)  P.  W. Forsbergh,  Jr,,  Phys.Rev._76,  1 187(1949). 

4)  R.G. Rhodes,  Acta  Cryst.  _2,  417(1949). 

5)  Terminology  according  to  the  "Internationale  Tabeilen  zur  Bestimmung  von 
Kristalstrukturen,"  Berlin,  Gebruder  Borntraeger,  1935. 
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1 * o 

11.  Tetragonal}  space  group  C^v  = P4mm|  aQ  = 3.9860,  cq  = 4.0259A 

Cell  contents}  BaTiO^ 

Ba  in  (a),  0,  0,  0 (Origin) 

Ti  in  (b),  1/2,  1/2,  (1/2  + zTi) 

Olin(b),  1/2,  1/2,  £q j 

*2  Oil  in  (c), -0,  1/2,  (l/2  + zqTi)}  l/2,  0,  (l/2  + *0II) 

Three  parameters}  stable  between  120°  and— 5°C. 

1 A Q 

III.  Orthorhombic;  space  group  C^v  = Ammj  at  4 C,  aQ  = 4.0185, 

b = 3.9860,  c = 4.0162^ 
o o 

Cell  contents:  2BaTiO^ 

2  Ba  in  (a)  0,  0,  0 (Origin);  0,  1/2,  l/2 
2 Ti  in  (b),  l/2,  0,(l/2+zTi)}  l/2,  l/2,  zTi 

2 OI  in  (a),  0,  0,(l/2  + Zq^);  0,  l/2,  Zq^ 

4 Oil  in  (e),  l/2,(l/4  + yOII),  (l/4  + zon);  1/2,  (3/4  + *OII)l 
l/2,(3/4-yon),  (1/4  + zon);  1/2,  ( 1/4  - yQ1I), 

(3/4  + Z0II) 

Four  parameters;  stable  between —5°  and— 90°C. 

IV.  Rhombohedral;  space  group  C^v  = R3m{  aQ  = 4.035^; 

a = 90°  19* 

Hexagonal  cell  contents;  3BaTiO^ 

3 Ba  in  (a),  0,  0,  0 (Origin);  1/3,  2/3,  1/3;  2/3,  1/3,  2/3 
3 Ti  in  (a),  0,  0,  (l/2  + z T-);  l/3,  2/3,  (5/6  + zTi); 

2/3,  1/3,  (1/6  + zT.) 

9 O in  (b),  (5/6  +xq),  (1/6-xq),  zq;  (1/6  + xq),  (5/6  - xj, 
(l/3  + zq){  (1/2  + xQ),  (1/3  - xo),  (2/3  + zQ){  (5/6  + xq), 
(2/3  + xQ),  zQ;  (1/6  + xQ),  (1/3  + xQ)(  (l/3  + zq); 

(!/2  + xo),  xq,  (2/3  + zQ);  (1/3 -xq),  (1/6-xq),  zo; 
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(2/3  - xq),  (5/b  — xq),  (l/3+zo)j  -xq,  (1/3  - xq), 

(2/3  + zQ) 

Three  parameters}  stable  below— 90°C. 

All  these  structures  are  based  on  the  familiar  perovskite  type,  with 
the  pseudocube  edge  always  close  to  4,£.  The  exact  variations  of  cell  dimen- 
sions of  all  the  phases  have  been  accurately  measured  by  Megaw^  and  Rhodes, ^ 
The  attempt  to  determine  exactly  the  three  parameters  of  the  tetragonal  phase 
II  (at  room  temperature)  is  described  in  the  following  sections. 


3.  The  Tetragonal  Structure  Factors 

While  the  parameters  z^,  Zqj  and  a*  described  above  for  the 

o 

tetragonal  phase  are  small,  they  are  appreciable  (of  the  order  of  0.1  A),  and 
one  would  believe  at  first  that  the  ordinarily  sensitiv  /stal- structure  methods 
would  measure  them  in  a straight-forward  manner.  The  intensities  of  the  dif- 
fracted beams,  which  constitute  our  ultimate  experimental  data,  are  a func- 
tion of  the  square  of  the  structure  factors  F which  express  the  manner  in 
which  X rays  scattered  by  the  individual  atoms  are  combined  into  the  beams 
scattered  by  the  crystal  as  a whole. 


(hk/ ) 


= 2 f + kyj  + 

j J 


= 2 cos  2Tr(hx,  + ky.  + £z.)  + i £ f,  sin  2ir(hx.  + ky.  + i z.)  (1) 
jJ  J J J j J JJ  J 

= aj + 1 f fj  y • 

in  which  h,  k and  / are  the  order  numbers  of  the  diffracting  plane  in  the  three- 


dimensional  grating  of  the  crystal;  Xj,  y^  and  z,  are  the  co-ordinates  of  the  j 
atom  in  fractions  of  the  unit-cell  edges;  and  f.  is  the  scattering  power  of  the  j 

J 

6) H.D.Megaw,  Proc. Phys. Soc,  5*8,  133(1946). 

7)  R.G. Rhodes,  Acta  Cryst.  4,  105(1951). 


jth 


th 


atom  for  the  particular  reflection  (hkl)  (a  function  of  the  Bragg  angle  2 6). 
The  summation  is  taken  over  all  the  j atoms  in  the  unit  cell.  The  cubic  unit 
cell  above  1 20°C  contains  five  atoms  with  the  co-ordinates: 


j 

Ba 

i 

0 

a* 

£ 

Ti 

1/2 

1/2 

1/2 

O 

1/2 

1/2 

0 

O 

0 

1/2 

1/2 

o 

1/2 

0 

1/2 

If  these  co-ordinates  are  substituted  in  Eq.  (1),  the  b.  terms  vanish 

J 

and  F is  determined  as  shown  in  Table  1,  according  as  h and  £ are  odd  or 
even,  respectively. 


Table  1 

Structure  Factors  for  (hO t ) Reflections  For  Cubic  Barium  Titanate 


h odd 

h even 

£ odd 
£ even 

fBa  + fTi  “ fO 
fBa  ~ fTi  “ fO 

fBa  " fTi  ” fO 
fBa  + fTi  + 3fO 

We  now  consider  only  reflections  with  k = 0,  since  these  are  the  only  ones 
measurable  in  this  study  (cf.  Section  4). 

In  the  tetragonal  crystal,  the  unit  cell  contains  five  atoms  with  the  fol- 
lowing co-ordinates: 


j 

X. 

_2 

h 

z. 

J 

Ba 

0 

0 

0 

Ti 

1/2 

1/2 

1/2  + zTi 

01 

1/2 

1/2 

Z0I 

Oil 

0 

1/2 

\fl  + zon 

on 

1/2 

0 

i / 2 + zon 
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The  new  Zj  co-ordinates  modify  the  structure  factors  in  Table  1 according  to 
Eq.(l),  by  multiplying  the  fj  terms  by  cos  2irizj  and  adding  an  imaginary  com- 
ponent in  fj  sin  2ir/zj.  But  since  the  co-ordinate  parameters  zTi,  zQI  and 
*011  ***  near  aer°i  is  apparent  that  the  absolute  values  of  the  structure 
factors,  | P'|,  will  change  very  little  for  relatively  large  changes  in  the  para- 
meters. Thus,  the  structure  factors  of  Table  1 are  a very  good  approximation 
to  those  of  the  tetragonal  structure.  In  fact,  the  approximation  is  so  good 
that  the  parameters  Zj  cannot  be  determined  by  ordinary  means. 


4.  Range  of  X-Ray  Measurements 

In  order  to  get  the  most  out  of  the  function  of  Eq.(l),  the  value  of  t in 
the  Bragg  reflections  measured  must  be  extended  as  far  as  possible.  In  ad- 
dition, it  is  desirable  to  have  as  many  Bragg  reflections  measured  as  possible 
so  that  random  errors  will  be  minimized.  To  this  end,  MoKa  radiation  was 
used,  giving  a range  of  (sin  0)/X.  from  0 to  1.41.  The  crystal  is  sufficiently 
hard  so  that  reflections  are  not  markedly  damped  by  temperature  vibration, 
and  are  still  appreciably  strong  even  at  the  limits  of  this  range.  The  number 
of  (hk i)  reflections  in  this  range  is  about  6,000,  but  because  of  the  difficulties 
of  correcting  for  absorption,  we  must  be  content  with  measuring  only  the  (h0^) 
reflections,  of  which  there  are  405.  With  the  Geiger  counter  instrumentation 
we  have  been  able  to  measure  356  of  these,  and  these  data,  after  averaging, 
yield  99  symmetrically  nonequivalent  Bragg  reflections. 


5 . Experimental  Source  of  the  Data 
The  essential  data  (99  F terms)  were  taken  from  a single  untwinned 
crystal  of  barium  titanate  at  room  temperature.  This  crystal,  whose  dimen- 
sions are  0.071  x0.160x0.193  mm,  was  taken  from  a preparation  made  by 
B . M.  Matthias.*^  The  crystals  were  grown  in  a flux  of  barium  chloride,  with 


8)  B.  Matthias,  Phys.Rev,  73,  808(1948). 
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a sm^dl  excess  of  barium  oxide  over  titanium  dioxide,  in  the  presence  of  one' 
mole  ratio  of  boron  oxide  to  barium  titanate  to  minimize  twinning.  The  crystal 
was  mounted  on  a glass  fiber  and  placed  on  the  standard  Weissenberg  camera 
for  X-ray  study,  A Geiger  counter  was  used  to  measure  356  integrated  re- 
flections as  described  briefly  in  Section  6,  These  measurements  were  cor- 
rected for  Lorentz  and  polarization  effects  in  the  standard  manner.  Correc- 
tions were  also  made  for  absorption  and  secondary  extinction  of  the  X rays  in 
the  crystal  according  to  the  method  described  in  Section  7.  These  corrections 
may  be  summarized  in  the  following  expressions 


I = AVQ1 
o 

Q = KLp|F|2 


(2) 


where  I is  the  integrated  intensity j 1Q  is  the  intensity  of  the  incident  beam;  K 
is  a constant;  L is  the  Lorentz  factor,  which  is  a measure  of  the  time  in  which 
a set  of  planes  has  an  opportunity  to  diffract  as  the  crystal  rotates;  p is  the 
polarization  factor,  introduced  because  an  unpolarized  beam  is  partially  po- 
larized in  the  diffraction  process  to  a degree  dependent  on  the  Bragg  angle;  V 
is  the  volume  of  the  crystal;  A is  the  absorption  factor,  dependent  on  the  linear 
absorption  coefficient  n (a  function  of  I in  this  case,  where  secondary  extinction 
is  important,  and  different  for  each  reflection),  the  shape  of  the  crystal  and 
the  geometry  of  diffraction;  and  F is  the  structure  amplitude  from  which  it  is 
hoped  to  determine  the  details  of  the  structure  itself  through  Eq.(l). 

The  99  symmetrically  independent  structure  amplitudes  obtained  by 

averaging  the  356  measured  reflections,  together  with  Lp  and  A corrections 

2 2 2 

and  the  Bragg  parameter  s = (sin  9)  /\  , are  listed  in  Table  2.  The  final 
set  of  calculated  values  of  F as  derived  in  the  following  sections  is  also 
given. 
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Table  2 

X-Ray  Diffraction  Data  For  Tetragonal  Barium  Tetanate 


0.564 

0.771 

1.002 

1.268 

1.568 

0.0154 

0.0310 

0.0778 

0.157 

0.266 

0.407 

0.581 

0.806 

1.014 

1.281 

1.583 

0.0615 

0.0767 

0.124 

0.203 

0.312 

0.450 

0.626 

0.826 

1.065 

1.332 

1.628 

0.138 

0.153 

0.200 

0.280 

0.387 

0.527 

0.701 

0.903 

1.136 

1.407 

1.698 


!(**!•)  I/Lp 


|F|  A B 

calc.  AF  t calc.  calc. 


880 

3410 

1042 

1820 

403 

642 

118 

203 

49 

98 

32 


1.207 

0.178 

0.256 

0.422 

0.632 

0.891 

1.203 

1.557 

1.877 

1.997 

1.750 

1.172 

0.369 

0.418 


0.743 

0.994 

1.308 

1.640 

1.925 

1.982 

1.672 

1.071 

0.586 

0.626 

0.740 

0.920 

1.163 

1.459 

1.767 

1.979 

1,929 

1.534 

0.900 


253 
405  0.058 
215  0.212 
248  0.177 

203  0.234 
215  0.216 

199  0.232 

204  0.217 
198  0.208 
202  0.199 
198  0.180 

407  0.035 
215  0.144 
294  0.104 

205  0.202 
223  0.202 

200  0.219 

207  0.212 
198  0.218 

202  0.204 
198  0.197 
205  0.176 

208  0.126 
248  0.085 
205  0.148 
219  0.164 
200  0.194 
209  0.183 
198  0.217 

203  0.196 
198  0.188 
202  0.181 
198  0.179 


104.0 

24.2 
65.5 
20.9 

40.3 
.15.6 
27.8 

10.7 

18.7 

8.1 

12.8 


0 

5.4  0 
5.4  0 
3i6  0 
4.0  0 


0.8  -2.0 

.2  1.3 

.3  -1.7 


0 

-1.1 
.5 
.4 
.5 
1.0 
0.3  -2.3 
0.9  0.9 

■1.5  -9.6 
3.2  4.8 
0.9 
0.7 


104.0 

24.2 

65.5 
20.9 

40.3 

15.6 
27.8 


‘T .'zCi  -KT  v , 


2 

0 

3 

0 

4 

3 

5 

3 

6 

0i 

0.386 

0.500 

0.634 

0.810 

1.008 

1.250 

1.513 

0.384 

0.398 

0.448 

0.521 

0.629 

0.774 

0.951 

1.158 

1.385 

1.646 

0.546 

0.564 

0.610 

0.692 

0.797 

0.945 

1.121 

1.311 

1.548 

0.755 

0.771 

0.815 

0.897 

1.004 

1.141 

1.306 

1.513 

1.750 

0.982 

0.998 

1.042 

1.124 

1.232 

1.374 

1.543 

1.750 


0.844 

237 

0.879 

203 

0.985 

223 

1.160 

200 

1.390 

210 

1.661 

199 

1.905 

204 

1 .920 

198 

1.803 

202 

1.314 

198 

1.152 

200 

1.187 

215 

1.287 

200 

1.448 

209 

1.657 

199 

1.868 

204 

1.996 

198 

1.920 

202 

1.572 

198 

1.017 

206 

1.500 

208 

1.535 

199 

1.622 

207 

1.754 

198 

1.900 

204 

1.994 

198 

1.950 

202 

1.691 

198 

1.230 

203 

1.838 

198 

1.858 

204 

1.911 

1.973 

1.999 

1.927 

1.698 

1.308 

0.788 

2.000 

2.000 

1.990 

1.944 

•1.825 

1.596 

1.248 

0.796 


198  0.251 
203  0.230 
198  0.229 
202  0.220 
198  0.230 
203  0.232 
198  0.264 

203  0.275 
198  0.283  ' 
202  0.268 
198  0.269 

202  0.264 
198  0.261 

203  0.263 
198  0.284 


Table  2 — Continued 


ho£ 

2 

s 

l(rel.) 

I/Lp 

A 

If! 

obs. 
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009 

1.250 

61 

1.810 

198 

0.304 

8.6 

8.1 

0.5 

-21.0 

7.6 

-2.9 

109 

1.254 

160 

1.787 

202 

0.305 

13.8 

14.1 

-0.3 

11.5 

13.8 

2.8 

209 

1.304 

44 

1.715 

198 

0.298 

7.2 

8.0 

-0.8 

-19.7 

7.5 

-2.7 

309 

1.381 

150 

1.576 

202 

0.290 

12.8 

12.8 

0 

11.3 

12.6 

2.5 

409 

1.491 

56 

1.357 

198 

0.290 

7.3 

7.3 

0 

-19.2 

6.9 

-2.4 

509 

1.633. 

148 

1.055 

206 

0.287 

10.6 

11.2 

-0.6 

10.8 

11.0 

2.1 

0010 

1.538 

172 

1.274 

203 

0.315 

11.9 

12.4 

-0.5 

7.0 

12.3 

1.5 

1010 

1.545 

72 

1.239 

198 

0.321 

7.5 

7.5 

0 

-11.5 

7.3 

-1.5 

2010 

1.593 

171 

1.137 

205 

0.316 

11.2 

12.1 

-0.9 

6.7 

12.0 

1.4 

3010 

1.672 

73 

0.964 

198 

0.325 

6.7 

7.2 

-0.5 

-11.2 

7.1 

-1.4 

6.  Geiger  Counter  Intensity  Measurements 
Special  apparatus  and  techniques  were  devised  to  measure  diffraction 
intensities  given  by  the  crystal  mentioned  in  Section  5, with  the  greatest  pos- 
sible accuracy.  Much  time  was  expended  in  an  attempt  to  use  film  photo- 
metric methods  to  evaluate  intensities  but  these  were  finally  abandoned  be- 
cause it  was  felt  that  random  errors  in  the  measurements  were  intolerably 
high.  The  Geiger  counter  was  finally  resorted  to  as  the  most  sensitive  and 
accurate  instrument  for  measuring  diffraction  intensities. 

The  equipment  used  consisted  of  an  adaptation  of  the  standard  Norelco 
Wide-Range  X-Ray  spectrometer  and  circuits  to  a standard  Supper  Weissen- 
berg  camera.  The  overall  assembly  is  shown  in  Fig.  2.  The  Geiger  counter 
is  attached  to  a bracket  pivoted  on  an  axis  coincident  with  the  spindle  axis  of 
the  Weissenberg  camera  as  shown  in  Fig.  3.  The  angular  position  of  the  Geiger 
counter  is  indicated  with  a vernier  on  a large  graduated  circle.  The  tilting  ad- 
justment shown  in  the  illustration  is  intended  for  equi-inclination  settings  on 

upper  levels,  not  used  in  this  investigation.  A back-geared  cam  and  lever  ar- 

» 

rangement  fastened  to  the  spindle-dial  end  of  the  Weissenberg  camera  per- 
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mitted  an  accurately  uniform  oscillation  of  the  crystal  through  a 0°  to  3°  range 
(continuously  selectable),  and  the  operation  of  a switch  to  start  and  stop  the 
counting  circuits.  The  intensity  was  integrated  by  scanning  the  MoKa  peak  in 
the  spectrum  produced  by  the  diffracting  planes  under  study.  If  a 1°  oscil- 
lation of  the  crystal  was  used,  a Geiger  counter  slit  opening  corresponding  to 
ca.  3°  yvith  the  counter  at  a fixed  angular  position  was  ample  to  insure  that 
all  of  the  spectrum  produced  passed  into  the  tube. 

Typical  scans  as  recorded  on  the  strip  chart  are  shown  in  Fig.  4.  Count- 
ing was  started  and  stopped  at  a precise  moment  of  the  scanning  cycle,  and 
the  counts  were  accumulated  on  the  count  register  during  one  oscillation  cycle, 
with  one  scan  forward  and  one  backward  (1  to  2 degrees  oscillation,  256 
seconds  of  time).  The  total  of  counts  was  then  proportional  to  the  area  under 
the  strip  chart  record,  made  simultaneously.  The  effect  of  background  was 
measured  by  determining  the  background  height  from  the  chart,  and  sub- 
tracting this  value  (multiplied  by  a suitable  factor)  from  the  count  total.  The 
result  was  proportional  to  the  integrated  intensity  for  the  reflection  corres- 
ponding to  molybdenum  Ka  monochromatic  radiation.  This  background  cor- 
rection allowed  the  use  of  unfiltered  radiation  with  the  consequent  gain  in 
intensity.  An  attempt  to  use  crystal-monochromatized  radiation  was  found  to 
be  useless  because  only  16  out  of  99  reflections  could  be  detected  at  such  low 
intensity  levels. 

Considerable  care  was  taken  to  insure  the  accuracy  and  reproducibility 

* 

of  these  measurements.  It  was  found  that  close  regulation  of  the  primary  X- 
ray  beam  was  necessary}  an  X-ray  tube  space  current  regulator  was  used  to 
maintain  this  current  constant  without  drift  to  0.1  percent.  This  regulation 
together  with  0.1  percent  primary  voltage  regulation  was  found  to  be  adequate. 
A standard  reflection  was  remeasured  periodically  throughout  the  work.  In- 
tensity measurements  were  reproducible  to  about  1 percent  where  the  statis- 
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tical  error  would  be  about  0.3  percent.  Variations  among  equivalent  reflec- 
tions, resulting  probably  from  slight  variations  in  absorption  and  extinction 
conditions,  were  frequently  higher,  sometimes  as  much  as  5 percent,  and 
occasionally,  for  very  strong  reflections,  up  to  10  percent. 

The  integrated  intensities  measured  as  described  were  averaged  for 
the  99  nonequivalent  reflections,  and  are  listed  in  Table  2 under  the  heading 
I (rel).  These  numbers,  corrected  for  background,  are  the  count  totals  di- 
vided by  the  scale  factor  of  128.  Further  details  of  the  measurement  of  in- 

91 

tegrated  intensities  with  a Geiger  counter  are  described  elsewhere.  ' 


7.  Correction  of  Intensities  for  Absorption  and  Extinction 
Among  the  most  important  factors  which  influence  the  final  value  of  the 
diffracted  beam  intensity  are  absorption  and  extinction  effects  (A  in  Eq.  2). 
This  factor  is  a function  of  the  linear  absorption  coefficient  p,  the  shape  of 
the  crystal,  the  position  of  the  crystal  in  the  beam  and  the  Bragg  angle.  For 
this  crystal  it  was  derived  from  an  exact  evaluation  of  the  expression 


A'"* 

/ 


over  the  whole  volume  of  the  crystal,  where  t is  the  path  length  of  the  X rays 
into  and  out  of  the  crystal  for  the  volume  element  dv.  In  order  to  make  the 
calculation  feasible,  only  diffraction  beams  normal  to  the  prism  axis,  corres- 
ponding to  the  (hO l)  reflections,  were  considered.  Reducing  the  problem  in 
this  way  to  the  two  dimensions  of  the  right-prism  section,  it  was  possible  to 
subdivide  this  section  into  areas  for  which  the  integral  of  Eq.(3)  is  continu- 
ous, with  boundaries  determined  primarily  by  the  edges  of  the  crystal.  Two 
typical  subdivisions  are  shown  in  Fig.  5.  The  limits  of  the  exponential  integral 
in  Eq.  (3)  are  different  for  each  area  j,  so  that  Eq.  (3)  becomes 

9)  H.  T.  Evans,  Jr.,  Technical  Report  No.  39,  Philips  Laboratories,  Inc,, 

1951}  Rev.  Sci.  Instr.,  in  press. 


where  Z.  is  the  exponential  integral  of  each  subdivision  multiplied  by  p and 


s is  the  area  of  the  section  for  the  two-dimensional  case.  In  practice,  for 


each  reflection  the  integrations  Zj  were  set  up  and  carried  out  (4  to  9 per  re- 
flection), evaluated  and  summed.  These  tedious  calculations  were  carried 
out  for  99  reflections  using  the  calculated  value  of  p = 198  em"^. 


Fig.  5.  Typical  subdivisions  of  crystal  section  for  determination  of  X- 
ray  absorption  effects:  (a)  for  (203)  reflection;  (b)  for  (501)  reflection. 

Application  of  these  values  of  A to  the  measured  intensities  showed  by 
comparison  with  calculated  structure  factors  that  secondary  extinction  is  an 
important  factor.  According  to  this  effect  the  drain  of  the  energy  of  the 
primary  X-ray  beam  into  the  diffracted  beam  resulted  in  increasing  the  ap- 
parent absorption  coefficient  p by  an  amount  proportional  to  the  intensity  of 
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the  diffracted  beam.  The  effective  p becomes  under  these  conditions 

V ~ (5) 

where  Q is  a function  of  the  structure  factor  and  includes  the  Lorentz  and 
polarization  effects  (Eq.  2)  and  g is  a pure  number,  determined  experimentally 
for  a given  crystal.  An  attempt  was  made  to  determine  the  value  of  g by 
making  integrated  intensity  measurements  of  the  strongest  reflections  with 
the  wide-range  spectrometer,  using  a flat  powdered  specimen.  It  was  found 
that  preferred  orientation  effects  made  it  impossible  to  account  for  the  in- 
tensities obtained,  even  though  extinction  effects  are  presumably  absent,  The 
value  of  the  factor  was  therefore  determined  by  trial,  using  calculated  values 
of  Q.  Different  values  of  g were  used  to  determine  p for  12  strong  reflec- 
tions until  the  best  fit  between  calculated  and  observed  intensities  was  ob- 
tained. This  process  is  demonstrated  in  Table  3.  The  reliability  factor  R 


Table  3 

Determination  of  Secondary  Extinction  Jn  Barium  Titanate 
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(see  Section  8)  is  used  as  the  criterion  of  fit,  and,  on  this  basis,  the  value 
g = 5.06 x 10-5  was  chosen.  All  absorption  factors  A were  then  recalculated 
and  these  values  were  used  to  establish  the  final  observed  structure  factors. 
Table  3 lists  the  values  of  A so  obtained,  and  also,  under  p,,  the  effective 
absorption  coefficient  as  determined  by  the  secondary  extinction  effect,  which 
was  used  in  the  calculation  of  A. 

In  spite  of  the  rather  drastic  thermal  history  of  the  crystal,**)  the  use 
of  molybdenum  radiation,  and  the  excellent  overall  agreement  between  ob- 
served and  calculated  structure  factors,  certain  of  the  strongest  reflections 
are  anomalously  out  of  line  with  this  general  agreement.  For  example,  the 
(200)  and  (002)  reflections,  which  have  calculated  structure  factors  both  equal 
to  65.5,  are  observed  to  have  values  of  60.1  and  73.4,  respectively.  Such  an 
effect  must  result  fi  the  influence  of  an  appreciable  amount  of  primary  ex- 
tinction, This  effect  will  be  observed  in  regions  of  the  crystal  (within  the 
"mosaic  blocks"  of  the  normally  imperfect  crystal)  which  are  so  perfect  that 
the  doubly-reflected  Bragg  reflection,  180°  out  of  phase  with  the  incident 
beam,  will  tend  to  cancel  out  this  beam.  Since  there  is  no  way  to  account 
quantitatively  for  this  effect  under  the  present  conditions,  it  has  been  necessary 
to  exclude  the  seriously  affected  reflections  from  consideration  in  the  structure 
analysis.  To  keep  this  selection  otherwise  impartial,  all  reflections  (16)  with 
a value  of  (sin  0)2/\^  < 0.250  have  been  excluded,  thus  leaving  83  reflections 
to  constitute  the  basic  data. 

The  matter  of  extinction,  secondary  and  primary,  will  not  be  further 
treated  in  this  discussion.  These  phenomena  have  been  thoroughly  investi- 
gated by  ICdnzig*^)  with  results  which  have  a highly  significant  bearing  on 
textural  phenomena  which  occur  in  barium  titanate  single  crystals  as  they  are 
10)  W.Kanzig,  Helv.  Phys,  Acta  24,  175(1951). 
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heated  and  cooled  through  the  Curie  point  at  120°C. 

The  absorption  calculations  described  above  (and  in  more  detail  else- 
where* ^ can  be  carried  out  with  some  effort  in  cases  where  the  plane  of  dif- 
fraction is  normal  to  the  prism  axis  of  the  crystal,  but  become  prohibitively 
involved  for  inclined  reflections.  Therefore,  of  the  nearly  6,000  reflections 
of  the  type  (hk£)  available,  only  those  with  k = 0 can  be  corrected  for  absorp- 
tion. It  should  be  noted  also  that  the  crystal  must  be  treated  as  it  stands, 
since  any  attempt  to  change  its  shape  would  fill  the  crystal  with  glide  twin 
lamellae. 


8.  Calculation  of  Structure  Factors 


The  structure  amplitude  F for  reflection  (hki)  is  a function  of  the  atomic 
co-ordinates  Xj,  y^  and  and  the  scattering  factor  f^  of  each  atom,  according 
to  Eq.  (1).  The  scattering  factors  used  in  this  investigation  are  those  of 

5) 

Hartree,  ' and  Pauling  an 
temperature  factor,  thus 


51  51 

Hartree,  ' and  Pauling  and  Sherman  ' for  the  stationary  atom  modified  by  a 


fj  = f(T)j  * 


■B.s 

J 


(6) 


where  f^j  is  the  tabulated  scattering  factor  for  the  atom  j at  rest,  Bj  is  a 
constant  proportional  to  the  square  of  the  mean  amplitude  of  vibration  u of 
the  atom,  and  s = (sin  9)/X.  According  to  the  Debye -Waller  theory  this  re- 
lation is  approximately  B = 8tr2  u^.  The  complete  form  of  Eq.  (1)  will  contain 
12  free  parameters:  3 structure  and  9 anisotropic  temperature  parameters. 

The  evaluation  of  the  parameters  has  been  carried  out  by  means  of  a 
combination  of  trial-and-error  operations  and  least-squares  analysis.  The 
quality  of  the  given  structure  is  tested  by  comparing  the  observed  and  calcu- 
lated structure  amplitudes.  The  usual  criterion  used  in  statistics  is  the  sum 
of  the  squares  of  the  residuals,  or  quadratic  mean  error  c defined  as 

11)  H.  T . Evans,  Jr .,  Technical  Report  No .45,  Laboratory  for  Insulation  Re- 
search, 1951;  J.  Appl.  Phys.  23,  663(1952). 
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where  n is  the  number  of  observations  (reflections)  and  p the  number  of 
degrees  of  freedom  (free  parameters).  The  scale  of  | F| 0^s  Is  established 
by  multiplying  by  the  factor  Kl 

k=  s;  i*-t|bs 

It  is  the  quantity  ( JF*!  0^s  — l^*lcaic)^  which  is  minimised  analytically  in  least  - 
squares  analysis. 

Crystallographers  have  been  in  the  habit  of  using  another  criterion, 
the  reliability  factor  R,  defined  as  follows! 


R=  S||Flob<,-|Ficalc|-^  StF|oba  (8) 

where  the  scale  for  |F|obfl  is  established  by 

k=  z inttlCT  s ui*.  • 

The  factor  R is  dimensionless  and  therefore  useful  to  some  extent  in  com- 
paring the  relative  degree  of  refinement  of  structures  of  different  crystals 
and  substances.  As  a general  rule,  a value  of  R as  low  as  0.20  is  considered 
satisfactory  proof  that  a given  structure  is  nearly  correct.  In  recent  years 
structure  determinations  giving  reliability  factors  of  0.12  and  0.15  have  be- 
come more  common  and  are  considered  extremely  good.  Only  by  means  of 

ionization  chambers  and  Geiger  counters  has  it  been  possible  in  a few  cases 

1 2) 

to  obtain  values  of  R as  low  as  0.06  (see,  for  example,  Cochran,  and 

13) 

Brill,  Grimm,  Hermann,  and  Peters.  ' 

1 2)  W.  Cochran,  Acta  Cryst.  _4,  81(1951). 

13)  R.  Brill,  H.G.  Grimm,  C.  Hermann  and  C . Peter s,  Ann.  Physik  ^1* 
393(1939). 
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9.  Factors  Influencing  the  Values  of  the  Structure  Amplitude 
The  values  of  the  observed  structure  amplitude  are  affected  by:  (1) 

The  relative  location  of  the  atoms;  (2)  the  scattering  factor  of  each  atom; 

(3)  errors  in  the  experimental  determination  of  the  F values.  We  try  to  re- 
duce the  latter  factor  as  much  as  possible  in  order  to  be  able  to  learn  more 
about  the  first  two.  ^ _ 

The  scattering  factor  fj  is,  to  a first  approximation,  a function  only 

of  the  Bragg  parameter  s (i.e.,  spherically  symmetrical),  and  is  similar 

•> 

-ks“ 

in  form  to  the  function  Ze  (Z  = atomic  number).  The  Gaussian  function 
is  modified,  of  course,  by  the  modern  statistical  mechanical  treatment  and 
has  been  calculated  by  Pauling  and  Sherman,  Hartree  and  others  for  all  atoms 
in  free  space  at  rest.  In  the  environment  of  the  crystal,  the  scattering  fac- 
tor will  be  altered;  first,  by  the  polarisation  of  the  atom  due  to  the  electric 
field  in  which  it  exists  and  bond  formation  with  other  atoms;  and  second, 
by  the  heat  motion  of  the  atom;  both  these  effects  will  be  aspherical.  The 
temperature  effect  makes  the  electron  cloud  of  the  atom  generally  more  dif- 
fuse, and  consequently  lowers  the  scattering  factor  at  high  Bragg  angles 
approximately  according  to  Eq.(6).  The  polarisation  of  the  atom,  on  the  other 
hand,  will  affect  almost  exclusively  the  outer  shell  of  electrons,  and  there- 
fore will  alter  the  scattering  factor  primarily  at  low  Bragg  angles,  while 
these  electrons  make  almost  no  contribution  to  the  high-angle  scattering. 

We  may  note  in  passing,  that  one  may  expect  to  obtain  a higher  degree  of 
agreement  between  observed  and  calculated  structure  factors  in  structures 
containing  heavy  atoms,  where  a higher  proportion  of  the  X-ray  scattering 
arises  from  dense,  closed-shell  atom  cores,  than  in  the  case  of  organic 
structures,  say,  where  a higher  proportion  of  electrons  are  influenced  by 
polarization,  since  polarization  is  the  hardest  factor  to  define  in  terms  of 
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the  structure  factor. 

In  the  case  of  barium  titanate,  because  of  the  considerable  influence 

of  primary  extinction,  it  was  necessary  to  discard  from  consideration  all 

2 2 

reflections  for  which  (sin  3)  /\  < 0.250.  This  means  that  with  83  data 
from  which  to  determine  12  parameters,  we  may  expect  to  learn  in  some  de- 
tail the  structure  (atomic  displacements)  and  temperature  motions  of  barium 
titanate,  but  will  not  be  able  to  learn  anything  directly  about  the  electronic 
polarizations  of  the  atoms  or  the  distribution  of  electrons  between  the  atoms. 

10.  First  Approximate  Model 

The  first  step  in  the  analysis  of  the  83  data  was  to  make  the  usual  as- 
sumption that  all  atoms  vibrate  spherically  wkh  the  same  amplitude,  giving 
a structure  defined  by  1 temperature  parameter  B,  and  3 structure  par- 
ameters zTi,  Zqj,  Zqjj*  By  trial,  each  was  varied  in  turn  by  small  inter- 
vals, a set  of  structure  factors  calculated  at  each  step  and  compared  with 

the  observed  set  by  the  criterion  of  the  reliability  factor  R (Eq.8),  With 

"16  2 

zTi  = zOI  = zOU  = ® (cu^ic  structure)  and  B = 0.30x10  cm  , the  reli- 
ability factor  was  determined  as  0.089  2,  already  a very  low  value  as  ordi- 
nary crystal  structure  standards  go.  After  28  trials,  the  lowest  value  for 
R for  this  model  was  found  for  the  following  parameters:  zT.=  0.015  cycles 
= 0.061.&,  Zqj  = — 0,024  cycles  = — 0.081A,  - — 0.020  cycles  a — 0.097A, 

B = 0.30£2,  R = 0.0473.  This  value  for  R was  considered  to  be  so  unusually 
low  that  the  results  were  published  in  a note*^  before  further  refinement  of 
the  model  was  carried  out.  It  may  be  noted  that  the  value  of  B is  quite  low, 
corresponding  to  a hard,  covalently  bound  crystal  with  a high  melting  point, 
in  contrast  to  molecular  crystals  (such  as  organic  crystals),  where  B com- 
monly may  be  as  high  as  2.5&2. 

14)  H.  T.  Evans,  Jr Acta  Cryst.  4,  377(1951). 
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11,  Anisotropic  Temperature  Motions  in  Barium  Titanate 
The  next  step  in  the  refinement  of  the  model  was  the  assumption  of 
different  temperature  vibration  amplitudes  (still  spherically  symmetrical)  j 

I 

of  each  type  of  atom,  thus  increasing  the  number  of  temperature  parameters  * 

to  four  j Ba&,  BTi,  Bor  BOir  Startin8  with  the  previous  model  in  which  j 

all  these  parameters  are  taken  as  0.  30i2.  appropriate  corrections  ABgg)  ( 

Ab^,,  , ABqj,  and  A B^  were  calculated  by  the  method  of  least  squares 
(Section  12),  to  give  a 7-parameter  structure  for  which  R = 0.0455,  a slight 
improvement  over  the  first  structure.  The  final  model  adopted  assumes 
anisotropic  temperature  motion  in  all  atoms,  and  is  defined  by  3 structure 
and  9 temperature  parameters.  Assuming  that  Bj  varies  elliptically  with 
direction,  the  temperature  parameters  are  designated  as  Qga»  PBa» 

PTi*  °oi'  Pol’  for  the  three  atoms  on  the  fourfold  symmetry  axes,  where  a 
corresponds  to  the  [a] -direction  vibration  and  (S  to  the  [c] -direction  vibra- 
tions; and  Poii»  Y0n  for  the  oxygen  atom  on  the  two-fold  rotation  axis, 

where  Qqjj  corresponds  to  the  [a] -direction  vibration  normal  to  the  titanium- 

oxygen  bond,  Yqjj  to  the  [a]  -direction  vibration  parallel  to  the  bond,  and  ^OII  I 

to  the  [c] -direction  vibration.  These  parameters,  which  are  shown  in  Fig.  6, 

appear  in  the  temperature  factor  in  place  of  B,  in  Eq.  (6),  as,  for  example, 

J 

-(aBa  c°s2$  + PBa  sin2$)s2 

fBa  = f(T)Ba  e 

where  4>  is  the  angle  between  an  (h0$)  diffracting  plane  and  the  a crystal- 
lographic axi3.  In  setting  up  the  structure  factors  it  is  now  necessary  to 
consider  the  two  Oil  atoms  separately,  since  their  temperature  motions  con-  ' 

t 

t 

i 

tribute  in  a different  way  to  the  amplitude.  The  complete  expression  for  the 
structure  factor  is  now  written 

I 

i 

i 

i 

i 
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.(.Bl  CO.2*  + fBa 


£(T)Ba 


-i.Tl  =«2+ + ^Ti  si°2*)s2  f(I)Tl[(i;)cos  ut  *«] 


1 

«»2+ + poi  ■ta**»*  «jT)0  R;i)«* *•*»«] 
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l odd  - + 
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function  (or  the  , tractate 


amplitude,  T,  in  term,  o£ 


We  have  defined  a 

l2par,met'rEt  a Pb  , cTi^Ti-or  »or  *«*■  "«r  W’ 

r = t(*Ti'  *or  Z0U'  ».• 

Approximate  valae,  oi  each  of  the. e par 
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Eq.(9),  can  be  calculated  a set  of  | F | ca^c  values  to  be  compared  with  the  ob- 
served set  of  data  |F|  ^ can  k*  shown  analytically  that  these  parameters 

may  be  altered  by  correction  terms  Azqj,  Az^,  Aafia,  • • • « etc., 

. -l  . ' - . 

such  that  the  standard  deviation  will  be  a minimum  when  the  correction  terms 
are  derived  from  a series  of  simultaneous  equations  of  the  formi 

^ A»Ti  + |^J  A‘0I  + 14  + ' ' ‘ ' * ■ lF'^  = A'F1, 

We  may  set  up  as  many  of  these  so-called  "normal"  equations  as  there  are  par- 
ameter corrections  to  determine.  All  of  the  data  can  be  used  by  forming  co- 
efficients for  the  correction  terms  by  summing  cross-products.  Thus,  if 


•airi  . . aiFi  _n  ain  _ 


= C,  etc., 


then  the  normal  equations  may  be  formed  as  follows: 

(SA2)AzTi  + (Sab)  Azqi  + (Sac)  AzQn  + . . . * (Sax  A|f|) 

(Sba)  Azt1  + (Sb2)  Azoi  + (2bc)  Az0II+ . . . =(2bxA|f|)  (10) 

(2ca)Azt.  + (2cB)AznT  + (Zc2)  AznTT  + . . . = (2cxA|f|)  . 


The  solution  of  these  symmetrical  equations  should  lead  directly  to  the  set  of 

parameters  which  will  give  the  best  fit  with  the  observed  data  since  our  ap- 

1 5) 

proximate  model  is  presumably  already  very  near  to  the  best  one.  ' 

The  differentials  used  in  the  least-squares  treatment  of  barium  ti- 
tanate  were  derived  as  follows: 


15)  For  further  treatment  of  least-squares  methods  see  A. C,  Worthing  and 
J.Geffner,  "Treatment  of  Experimental  Data,"  New  York,  John  Wiley  and 
Sons,  Inc.,  1943}  E . T . Whittaker  and  G.  Robinson,  "Calculus  of  Observa- 
tions," London,  Blackie  and  Son,  Ltd.,  1944. 


F * A + i B 
lF|2  - A2  + B2 

2 ir|';2j£L  -2a.|A  *?b-|S 

a|Fl  . A 3A  B 3B 

d*  |f|  If! 


cos  + sin  \J/ 


3B 

a* 


> 


f 


(ii) 


where  \ft  is  the  phase  angle  of  the  noncentrosymmetric  structure  factor. 

In  the  next  stage  of  the  analysis  of  barium  titanate,  the  13  differentials 
were  calculated  for  the  83  (hOl)  reflections,  their  cross-products  summed  to 
form  the  normal  equations  according  to  Eq.(10),  end  these  equations  solved 
by  the  method  of  Crout.^  Until  the  last  stages  of  the  work,  all  these  oper- 
ations were  accomplished  with  a hand  computing  machine,  and  were  checked 
by  a complete,  independent  recalculation.  The  results  of  the  first  complete 
calculations  proved  to  be  inconclusive  for  reasons  to  be  described  in  the 
next  section. 


13.  Interactions  between  Parameters 
When  an  attempt  was  made  to  solve  the  13  normal  equations,  results 
were  obtained  which  were  obviously  without  any  significance,  i.e.,  the  cor- 
rected parameters  gave  a standard  deviation  much  higher  than  the  previous 
set.  A study  of  the  differential  coefficients  shows  that,  for  the  values  of 
the  parameters  used,  certain  pairs  of  the  derivatives  are  nearly  equal  in 
value  for  all  |F|  terms,  except  for  a common  factor.  This  means  that  all 
the  values  change  in  the  same  way  with  changes  in  either  one  or  the  other  of 
the  pair  of  parameters.  There  are  4 such  pairs  of  derivatives  among  the  13: 
d|Fi/3aTi  and  3|F|/3pTi?  3 |F|/3  zQ1  and  3 Jf|/ 3P0I>  3|F|/ d*OII  and 
3|F|/dpQjjj  3|F|/3a^,.  and  3|F1/3qqjj.  Thus,  the  atomic  displacement 


16)  P.Crout,  Trans.  Am.  Inst.  Elect,  Eng.  60,(1941). 


parameter  z and  the  [c^  -direction  temperature  parameter  for  the  atom  p, 
have  an  effect  in  the  function  Eq.(9)  which  does  not  permit  their  separation. 

In  addition,  the  [[a] -direction  vibration  parameter  for  titanium  cannot  be 
distinguished  from  that  of  oxygen  II  which  is  superposed  on  it  in  the  (hQ t) 
projection. 

The  interaction  can  be  demonstrated  analytically  in  the  following  way. 

In  the  least-squares  analysis,  we  may  expect  to  have  difficulty  in  determin- 
ing correction  terms  for  two  parameters  $ if  the  function  varies  in  the  same 
way  for  each.  In  the  present  case  this  circumstance  will  obtain,  for  example, 
if  the  ratio  d |F|/3  *Tir  31*1/ d?Ti  is  constant  or  nearly  so.  It  can  be  shown 
that  for  the  set  of  parameters  which  characterize  the  barium  titanate  struc- 
ture, the  derivative 


9i 
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Ti 


9lF| 

3 ^Ti 


_ BP 

'n 


tends  to  vanish.  The  first  derivatives  have  the  following  form: 


tllL.- 2*1, 


■(qTi  cos24>  + PTi  sin2$)s2 


3z 


Ti 


'(T)T1  [(-+)  **"  2,r<ITI  c*‘f 
+ (”  * ) cos  sin 


9IFI  2 . 2,  -(«Xi  + PTl  sin2<i.)s2 

-p—.-s  sin  *s  i(T)Tj 


(^.  + ) C0S  Zfl'/z.j,.  C0S  ^ 

+ ( — + ) sin  2*£zTi  sin 


2 6l 


2*  _ 


sin  4>  = 


4 sm  8 


, and  £ is  the  index,  then 


If  we  note  that 
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p . 8tt 
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( + “ y sin  2iri  Ejj  cos  \f/  +(  “^  ) cos  Zv£z^^  sin  \jf 
( + “ ) cos  2ir /*Ti  cos  \j/  + ( * ” ) sin  2tt£ zTi  sin  ^ J 


Sir 


ic#2  Y 

where  x is  the  numerator  and  y the  denominator  of  the  bracketed  fraction,  then 


3 P . 8ir 

ITT7 


l/l  3 x x dy  \ 1 * 

7ly  17"7--ff)-py 


When  Zj  is  small,  then  approximately 
cos  2wi*Ti  = 1 

COS  [f/  a 1 

3/3i  (cos  l// ) = 0 

x = (IT)  8in2ff  £ aTi  + (”  *)  Sin  ^ 


■(:;) 


i 


sin2trizTi  3 2ir^*Ti 


(2lrzTi  COfl  2nieTl  sin  \J/)-  ^ (ain2ir^zT.-  sin  t//) 


Since,  when  the  values  z.  are  all  small,  sin  ^ is  small  and  changes  slowly 

J 

with^,  the  derivative  dP/9£  is  seen  to  have  a negligible  value,  at  least 
for  titanium  and  oxygen.  When  BP/&£  approaches  zero,  the  determinant 
of  the  coefficients  of  the  normal  equations  tends  to-  vanish,  and  consequently 
the  error  in  the  separate  determination  of  the  variables  involved  in  P tends 


to  become  infinite. 
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Strictly  speaking,  if  we  do  not  arbitrarily  fix  the  origin  of  co-ordinates, 
we  should  introduce  a 13th  parameter  Zg  , in  which  case  we  will  find  that 
d|F|/  and  d|F|/  tend  to  interfere  with  one  another.  Thus  it 

appears  that  only  3 parameters  can  be  established  unambiguously}  ag&, 

aOI  and  Y0I1‘ 


The  least- squares  analysis  of  the  complete  model  (12  parameters)  next 
was  carried  out  with  8 parameters  assuming  that  one  of  each  of  the  inter- 
fering pairs  was  fixed.  Thus,  the  z parameters  obtained  from  the  original 
trial-and-error  study  were  assumed  to  be  correct  (and  also  ciqjj),  cor- 
rections were  determined  for  the  remaining  temperature  parameters.  The 
resulting  set  of  parameters  gave  a considerably  improved  agreement  between 
observed  and  calculated  structure  amplitudes,  with  R = 0.0372.  The  values 
obtained  for  the  temperature  vibrations  did  not  seem  at  all  reasonable  (see 
Table  4,  F 30),  however,  and  it  seemed  likely  that  the  assumption  the  z 
parameters  were  correct  must  be  abandoned.  The  calculation  of  the  independent 
parameters  gave  results  which  are  reasonable,  considering  only  the  binding 
forces  between  the  atoms.  Assuming  that  temperature  vibration  will  be  less 
in  the  direction  of  a bond  with  an  adjacent  atom  than  in  other  directions,  we 
expect  the  1 2-co-ordinated  barium  atom  to  be  quite  spherical  in  its  vibration, 
while  oxygen,  being  fairly  tightly  bound  to  titanium  in  one  direction,  should  be 
rather  anisotropic  in  its  vibration.  The  calculation  yields  the  following  results; 

GBa  = °*27 
PBa  = 0.28 

Gqj  = 0.90  (perpendicular  to  bond) 

Yqu  = 0.60  (parallel  to  bond). 
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Table  4 

Trial  Structures  Of  Tetragonal  Barium  Titanate 


F 

*Ti 

*OI 

m 

PBa 

aTi 

PTi 

°Sl 

Poi 

ES 

pon 

ESS 

R 

H 

KB 

■ 

am 

EH 

0.30 

.0892 

1.838 

20 

.015 

-.024 

-.020 

0.2 

\0 

.0473 

1.158 

29 

.015 

-.024 

-.020 

0.29 

0.39 

0.88 

. . 1 

0.49 

.0455 

1.053 

30 

.015 

-.024 

-.020 

.27 

.28 

.53 

.21 

.90 

.08 

.49 

.07 

.60 

.0375 

0.861 

42 

.012 

-.026 

0 

.27 

.28 

.46 

.30 

.90 

.50 

.90 

.90 

.60 

.0372 

0.924 

Kan- 

.014 

-.032 

■ 

0.48 

0.13 

0.13 

0.48 

.123 

We  make  the  following  adjustments] 

1.  Let  (Iqjj  = ^OII  = since  these  all  correspond  to  directions 
normal  to  O-Ti  bonds. 

2.  Let  = 0.50,  a little  smaller  than  because  the  01  Ti 

bond  is  presumably  somewhat  stronger  than  OIl-Ti. 

3.  Since  is  raised  to  0.90,  we  decrease  a^,.  to  0.46  in  ac- 

cordance with  the  relative  values  of  their  differentials  cal- 
culated from  Eq. (9).  -- 

4.  Let  = 0.30,  allowing  for  the  stronger  bond  in  this  direction. 

Then,  by  triad  and  error,  the  as  parameters  were  redetermined  to  give  a 
minimum  value  for  R which  was  substantially  equal  to  that  obtained  for  the 
previous  model.  This  required  rather  large  adjustments  of  the  structure, 
and  the  result  is  shown  under  F 42,  Table  4,  The  calculated  structure  ampli- 
tudes for  this  model  are  the  ones  listed  in  Table  3. 

Thus,  we  have  arrived  at  a final  result  which  is  definitely  ambiguous. 


t 


I 

/ 

s 


It  can  be  stated  simply  that  unless  data  can  be  obtained  extending  consider- 
ably beyond  the  present  limit  of  (sin  0)/\  (by  the  use  of  shorter  wavelengths), 
the  structural  displacements  in  barium  titanate  cannot  be  uniquely  determined. 
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The  'measurement  of  (hkl)  or  (hkO)  reflections  would  serve  to  determine  the 
parameters  aTi  and  separately,  but  the  z and  p parameters  will  inter- 
act for  all  reflections. 

Table  4 shows  the  course  of  the  structure  determinations  by  listing 
certain  trial  structures  (serial  number  under  F).  The  starred  parameters 
are  those  uniquely  determined.  Structure  42  represents,  we  suggest,  the 
most  reasonable  one  consistent  with  the  X-ray  measurements. 

15.  The  Work  of  KXnzig 

While  there  have  been  many  investigators  working  on  the  barium  ti- 
tanate  structure  problem,  only  one  has  presented  a report  containing  the 
kind  of  measurements  and  analysis  capable  of  solving  it.  The  work  of 
Kanzig1^  consists  of  the  measurement  of  X-ray  diffraction  intensities  from 
a large  single  crystal  plate  of  barium  titanate.  The  results  of  his  work  on 
the  structure  are  shown  in  Table  4,  There  are  several  important  aspects  of 
K*anzig's  report  which  must  be  considered  in  evaluating  these  results! 

1.  Only  18  reflections  were  measured:  8 (hOO)-type  reflections  which 

would  depend  on  o^j,  and  and  10  (ool)-type 

which  would  depend  on  zT.,  zQI,  zou,  PSa,  PTi,  pQI  and  pQn. 
These  measurements,  while  they  extend  to  the  same  range  of 
(sin  8)/\  as  those  made  here,  are  hardly  enough  to  lend  statistical 
weight  to  the  determination  of  the  parameters. 

2.  It  was  assumed,  quite  arbitrarily,  that  Zq^  = 0.  While  the  assump- 
tion is  apparently  nearly  correct,  no  study  was  made  to  establish 
this  result. 

3.  It  was  assumed  that 

QTi  = PTi  = aOI  = POI 


aBa  = PBa  = aOII  = POII  = ^OII  ’ 
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again  apparently  quite  arbitrarily.  Thus,  only  two  spherical 
temperature  parameters  were  considered. 

4.  No  numerical  data  of  any  kind  are  given  in  Kanzig's  report. 

It  would  appear  that  Kanzig's  structure  determination  is  quite  incomplete. 

On  the  other  hand,  as  mentioned  in  Section  7,  his  measurements  were  made 
ovier  a range  of  temperatures  emending  from  20°  to  350°C  and  indicate  some 
very  significant  diffraction  phenomena  in  the  vicinity  of  the  Curie  point. 

Here  all  extinction  effects  suddenly  disappear,  and  from  this  he  was 
able  to  deduce  very  useful  information  about  the  domain  structure  in  the  crys- 
tal, and  the  transformation  process  at  the  Curie  point.  \ 

\v  _ 

16.  The  Future  of  the  Barium  Titanate  Structure  Problem 

It  is  a considerable  disappointment  to  workers  in  the  theory  of  fer- 

roelectrics,  for  whom  any  new  Observational  data  is  of  the  utmost  importance, 

to  find  fliat  X-ray  diffraction  fails  in  this  case.  One  way  out  of  the  difficulty, 

as  has  been  mentioned,  is  to  use  shorter  wavelengths,  as  Grimsal17^  has 

tried  to  do,  but  even  if  sufficient  intensity  can  be  achieved  to  get  data  in  the 

high  ranges,  the  atomic  sea tteTing  functions  are  not  known  in  these  regions 

and  must  be  determined  or  calculated.  Otherwise,  recourse  may  be  made  to 

neutron  diffraction  where  the  scattering  factors  are  radically  different  in 

relative  magnitudes  and  do  not  damp  out  at  high  angles;  they  may,  therefore, 

give  a complete  picture  of  the  structure. 

Finally,  it  should  be  pointed  out  that  this  is  only  the  second  instance 

in  which  the  temperature  factors  have  interfered  with  the  X-ray  determination 

of  structure.  The  other  case  was  hydrocyanic  acid  HCN,  worked  out  by 

171  E.G.Grimsal,  "Structure  of  BaTiO^,"  U.S.  Atomic  Energy  Commission 
Report  No.  AECU-884,  Dec.  7,  1950. 
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18} 

Dulmage  and  Lipscomb,  ' which  came  to  the  writer's  attention  long  after  the 
situation  was  discovered  for  barium  titanate.  The  phenomenon  should  be 
watched  for  in  future  determinations,  especially  when  an  atom  is  slightly  dis 
placed  from  a symmetrical  position  in  a polar  structure. 
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